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1. The Motivation 



A common feature of the cuprate superconductors is that the coherence 
length is comparable with the lattice spacing [1]. Therefore the Cooper 
pairs are rather localized in the coordinate space and consequently can be 
regarded as boson degrees of freedom [2]. The superconductivity is thereby 
associated with the kinematical Bose-Einstein condensation and there must 
be uncondensed pairs in the normal phase. 

For a Bose-Einstein condensation to happen, the thermal wavelength of 

the bosons with mass mi at the transition temperature Tc, Xc — \j m ^ B Tc 
should be comparable with the inter-boson distance, /, so that quantum 
mechanical coherence takes place. Indeed, the ratio ^f- is 1.38 for an ideal 
Bose gas and 1.65 for Hell at the A-transition. For high Tc materials, this 
ratio can be extracted from the results of the fiSR experiment [3] of Uemura 
et. al. According to them, the transition temperature Tq of the under-doped 
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up to optimal-doped cuprates is inversely proportional to the square of the 
magnetic penetration depth at T = 0, Ao with a constant of proportionality 
universal for all cuprates, 

T C (K) = 0.25 x 10 5 4^ , (1.2) 
e A 

where m e is the electron mass in vacuum and the unit of m e /{e 2 \ t Q) is A -3 . 
Taking into account of possible contribution to the superfluid density from 
the fermionic component, we have 

1 4e 2 

" (L3) 

where c denotes the average distance between CuO<i layers and / stands for 
the average inter-boson distance within each layer. Combining (1.3) and 
(1.2), we end up with 

Y < 2.8, (1.4) 

which is consistent with the picture of Bose-Einstein condensation. 

With simply an ideal electron(hole) gas in chemical equilibrium with 
an ideal gas of Cooper pairs, we are able to explain Uemura's universal 
dependence of Tq on the superfluid density [4] . 



2. The Boson-Fermion Model 



The boson-fermion model was proposed independently in Ref. 5 and in 
Ref. 3. The grand Hamiltonian of the system reads 

h = j2 tp4s a ps + 

p,s p 

p,q 

where a^ s , at,, bp and &t are annihilation and creation operators of elec- 
trons (holes) and bosons, the subscript s denotes the spin orientation, and O 
is the total volume of the system. The momentum dependence of the cou- 
pling will determine the pairing symmetry, e. g., g = const — > s-wave 
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and g oc q 2 — q 2 , — > d-wave. To illustrate the main physics, we assume a 3-D 
isotropic jellium in which 

p 2 
2m/ 



and 



.2 



2m& 



w p = +ZT" + 2(i/ -a») (2-3) 



and a constant g, where mf(mb) is the effective mass of fermions(bosons), 
fx is the chemical potential of the system, 2u is the energy of a static boson 
relative to two static fermions (2u = uJp = Q — 2eJ0 = o), 
The conserved electric charge number is given by 

p,s p 

Because of the strong Coulomb repulsion, we assume that v > and the 
bosons are resonances. A dimensionless coupling can be defined as the ratio 
of the boson half width in vacuum, T/2, to the energy 2v, i.e. 



9 ~ <lv ~ 7T I 2 ) y/F {2 ' b) 



which serves as an effective expansion parameter. 

The physics of the model at weak coupling g « 1 and T = is deter- 
mined by the total number density of charges, 

n = nf + 2n& (2.6) 

in relation to the characteristic density 



(2^)3 , , 

v- 37r2 , 



where nf is the number density of fermions and n& that of bosons. n v cor- 
responds to a filled Fermi sea with Fermi energy v. If n < n v , the elec- 
trons (holes) will fill in the Fermi levels below v and there can only be vir- 
tual bosons through interaction condensating at the zero momentum bosonic 
level. The long range order in this case is of BCS type. If n > n v , the ex- 
cess electrons(holes) after filling up the Fermi level v prefer to combine into 
bosons which condensing at the zero momentum level and the long range 
order is of Bose-Einstein type. 
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The grand partition function of the model (2.1) is 

Q = Tre~P H (2.8) 
and the thermal average of an operator O is defined to be 

TtOe~P H , , 

< o >= -^qnr- (2-9) 



With a Bose condensate 



B = ^=< b p = >, (2.10) 



we may substitute 

b f = VtiBSpfl + Pp (2.11) 
into the Hamiltonian (2.1) and obtain H = Hq + Hi, with 

H = 20(z/ - /u)B 2 + ^ps a vs + u p4Pp 



p,s 



+gB J2(4\ a -?l + a -<fi a <ft) ( 2 - 12 ) 



q 



and H\ the rest of the terms of (2.1) which are of higher orders at weak 
coupling. B has been chosen real. The Hamiltonian (2.12) can be easily 
diagonalized through a Bogoliubov transformation and yields a fermionic 
spectrum 



E f =y/$+A* (2.13) 

with a gap energy A = gB. The thermodynamic potential within this ap- 
proximation reads 

In Q = 2/3(jjl - u)B 2 Q + H 1 + ^ PEp ) 



- ln(l - e -^ +2v - 2 ^). (2.14) 
v 

The thermodynamical equilibrium corresponds to the maximum of In Q at 
fixed (3 and fx, i.e. 

dlnQ\ 



dB 2 Jf3,n 



= (2.15) 



and 

( d 2 lnQ \ 
\(dB 2 ) 2 ) 

Combining (2.15) and the relation 



<0. (2.16) 



n = o(^rJ^' (2 - 17) 

we can solve for the gap energy A(T) and the transition temperature Tc in 
terms of the density n. 

At low density, n < n„, we find that 

A(0) = 8/j exp ( - 2 - " ~M (2.18) 



2 2 

with [i ~ ^ n ) 3 and that 

with 7 = 0.5772... the Eular constant. At higher density, n > n v , we obtain 
that 

k B Tc = 2^m b ( r ^^Y (2.20) 
V 2C(1) ' 

with C(3/2) = 2.612... and that 



A(T) = A(0)yi-(^) § (2.21) 

with A(0) = g-y/ (n — n v )/2. Therefore, the phenomenological model (2.1) 
provides a simple interpolation between a BCS condensation and a Bose- 
Einstein condensation. 



3. The Pseudo-gap 



A remarkable phenomenon of under-doped cuprates is the pseudo gap in 
their electron(hole) spectrum above T(j [6] . It has been found numerically in 
the boson- fermion model [7]. Here I shall present an analytical calculation 
[8] , which highlight the importance of the dimensionality of the system. 
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Starting with the retarded electron propagator 

l>O0 

Sr( P o,P) = -i J (to*** < {a p it), 4(0)} >, (3.1) 

where the time development of ap(t) and its conjugate follows from the total 
Haniltonian (2.1), the spectral function A(po,p) is given by its imaginary, 
i.e. 

A(P0,p) = -lmS R (p ,p), (3.2) 

The function A(po,p) gives rise to the probability density of a single electron 
(hole) excitation at energy po and momentum p and can be measured di- 
rectly by ARPES. For a BCS superconductor, A(po,p) is peaked at the quasi 

particle pole po = in the normal phase and at po = ±^Je"i+ A 2 below Tc- 

1). Two-dimensions: This case simulates the actual cuprate materials 
above Tc which consist of a set of parallel Cu02 layers. As long as T is 
not too close to Tq, the inter layer hopping can be neglected. With a jellium 
approximation, eg and remain given by (2.2) and (2.3) but with a 2D 
momentum p. On writing 

Sr(po,p) = =k, (3-3) 

the retarded self-energy function, T, R (pQ,p), is given to the one-loop order 
by 

v , - 2 f d 2 q N b (q) + N f (q-p) 

Z R (po,p)=g j7r-\2 7 — — (3-4) 

J (2tt) 2 po - u ? + z0+ 

with 
and 

e p + 1 

respectively. There is no Bose-Einstein condensation in a truly 2D system, 
the bosonic chemical potential 2{ji — u) = —5 vanishes at T = and the 
integral (3.4) diverges logarithmically. For an approximate 2D system and 
with a fixed carrier density, 5 << h^T within a considerable range of T « 
vjk^ above Tq [8], and the integral (3.4) is dominated at small q. This gives 
rise to a crude estimate that 

£*(po,£)~^- (3.7) 
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with 



A2 2r ~2 u rp, k BT 

A — — g vkbTYo. — — . 

7T 



(3.8) 



and r = nib/ntf, which leads to a perfect gap. A more refined calculation of 
the integration (3.4) was performed in [8] for po ~ k B T with the result at 
p = Pf (Fermi momentum): 



A(P0iP)\p=pf = - Im 



7r po - u(po) - iv(po) + u(S) ' 



(3.9) 



where 



HPo) = —9 v 

71 



k B Tsign(p - 8) 



In 



4r/j5 + (po - 8) 2 + \po - 8 



^Ar^8 + (po-5) 2 yjArjid + (p - 8) 2 - \p Q - 8\ 



lnr 



+ lf( e -P(Po-6)\ 
1 2 J 1 \i /•// 



and 



(3.10) 



v{po) = 9 v 

7T 



with 



nk B T 



^Ar^ + (po-5) 2 2 



' ixk B T 



f(z) = —/= I dx^fx Zt 

0r Jo (l + ze x ) 2 



(3.11) 
(3.12) 



and C(l/2) = -1.4604.... The function (3.9) is plotted in Fig. 1. 

The argument based on the logarithmic divergence can be generalized to 
higher orders in g and a Borel summation of the most divergent diagrams 
yields [8]: 



A(p Q ,p)\ P=PF = l -^e-^. 



(3.13) 



It consists of two peaks at po = ±A/V2 and a depletion of states at the 
Fermi level po = 0. 

2). Three dimensions: With three dimensional p and q, the simple argu- 
ment based on the logarithmic divergence does not work since the integral 
(3.4) converges as 8 — > 0. The evaluation of the integral for 3D case leads to 



iir r f°° 

^(po,p)\p= Pf = —g 2 k B T + -g 2 /i / dx 

4 7T Jo 



X 



e Pfi(x 2 -l) _|_ l 
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The spectral function at p 
2 (solid line) and D 



p F , with fiii = 10, r = 2, \nk B T/S = 8 and 
3 (dashed line), where x = Pa/k B T and y — 



A(p ,p)\ p=PF k B T. 



—g k B T 



In 



k B T 



5 + 



(po-s) 2 

4 M 



+ 21n(l + e 



(3.14) 



The corresponding A(po,p) function is also plotted in Fig. 1 for comparison 
and the pseudo-gap disappears. Since the 3D calculation is only to the one- 
loop order, this is not in contradiction with a weak pseudo-g inferred 
from the numerical solutions of truncated Dyson- Schwinger equations [9]. 

The dependence of the pseudo-gap can be tested by changing the sepa- 
ration between C11O2 layers, or by comparison with 3D strongly correlated 
superconductors, say fullerenes perhaps. 



4. DC Transport Coefficients 



On the normal phase DC transport coefficients (resistivity, Hall number 
and magnetoresistance), lies perhaps, the most severe disagreement between 
the observation and the conventional Fermi liquid theory. A large number 
of measurements converges to a universal temperature dependence above Tc 
for all cuprates [10]. The resistivity p depends linearly on T in contrast 
with the T 2 behavior of a good metal. The Hall number njj, instead of 
being a constant which equals to the actual carrier density in case of a good 
metal, increases linearly with T and reaches a value about a factor two of the 
chemically determined carrier density. The magnetoresistance, Ap/p being 
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proportional to T n (n = 3.5 — 4) is at entirely variance with the Kholer's 
law (which says that the ratio y/^r is temperature independent) for good 
metals. 

In the boson-fermion model, there are two kinds of charge carriers above 
Tq: electrons or holes (fermions) and uncondensed Cooper pairs (bosons). 
If the lifetime of the latter is longer than the relaxation time of the electric 
current, both of them contribute to the DC conductivity A Green's func- 
tion calculation is, however, very difficult since the simplifications based on 
incoherent superposition of quasi elastic collisions fail for bosons. Instead, 
we ask the following question [11]: Given a chemical equilibrium of fermions 
and bosons (g = in (2.1)) and a simple power dependence on T for the 
relaxation time of each kind of carriers, Is it possible to produce a reasonable 
fit of available experimental data? 

Consider the situation with an electric field E parallel to the CuOi plane 
and a magnetic field B = Bz perpendicular to it, the electric current density 
is given by [11] 

J a = a ab (B)E b (4.1) 
with a, b = x, y. The Taylor expansion in B of the conductivity tensor reads 

a ab (B) = (a + Aa)5 ab + a H e ab + 0(S 3 ), (4.2) 

where 

a=m + 4^y, (4.3) 
\ m,f m b / 

is the zero field conductivity, 



is the Hall conductivity and 

Act = 



r ' n/ +8^V5 2 , (4.5) 
mr / 



,3 



rrvj m b 



with rif and n b the number densities of each type of carriers, Tf and r b the 
corresponding relaxation times, and nif and m b the corresponding effective 
masses. We have assumed a parabolic band for bosons and a nonparabolic 
band for fermions. In case of the nonparabolic band, the effective mass nif 
is defined by 

1 1 r ^ [#(«£) (4-6) 



nif 2nf J (27r) 2 k 
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and the band correction factors r\ and rj' of (4.4) and (4.5) are given by 

/d?k 
j^VaViV jv b 5{ej:) (4.7) 



m 2 f 



2n f 
and 

m) 



/d k 
j^-^v a V b viViVjViS(eji>. (4.8) 



2n f 

It follows from (4.3)-(4.5) that the DC resistivity 



p=-J T -^+^y\ (4.9) 

e z V rrif / 



the Hall number 



2 (w +4 ^ 
°" D _ V m f m b J 



2 



= — B = ^- (4.10) 



2 



and the Kholer's ratio 



K =-p -2—^- = ^ H2+ I 7 ?— ) P • 

B z my L V my m^/ J 

(4.11) 

The approximate agreement with Luttinger theorem as observed by photo| 
emission ruled out any issue based on a large fraction of bosons and we expect 
that n b << nf. The temperature dependence of nj and n b can be neglected. 
The DC resistivity is dominated by fermions and we have 

r f ocT-\ (4.12) 

An explanation of this non -Fermi liquid behavior is suggested in [12]. If the 
factor rj is sufficiently small such that the denominator of (4.10) is dominated 
by the second term, a Hall number considerably larger than the actual carrier 
density emerges and the linear temperature dependence can be obtained if 
we assume 

nocT- 1 - 5 . (4.13) 

With both (4.12) and (4.13), the Kholer's ratio would be proportional to 
T~ 15 if the factor 7/ is also sufficiently small. For a nonparabolic band, it is 
possible to have small rj since it switches sign from a particle Fermi sea to a 
hole Fermi sea. But a small rj' may not be easy to tune since its integrand 
is positive definite. We are still not in the position to claim the triumph of 
the boson-fermion model in explaining the DC transport coefficients. 
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5. Two Theoretical Issues 



To justify the boson-fermion model as an adequate low density phe- 
nomenological model of certain strongly correlated electronic systems, the 
following issues have to be settled. 
5.1 Fermi distribution function 

Consider a purely electronic system on a lattice with electronic opera- 
tors ap s and ai g in the momentum representation. The Fermi distribution 
function is defined as 

s 

with | > the ground state of the system. With the Bogoliubov type of trial 
state [13] 

| >= exp [ d p( a -P,l a PA ~ a lj a -p,i) 1° >> ( 5 - 2 ) 
p 

we find that 

n = 2 sin 2 6p, (5.3) 
For a BCS ground state, we have 

cos 2%= —= €? (5.4) 



and 



sin 29 = ; (5.5) 

with A << ep and the corresponding is slightly smeared by A from that 
of an ideal Fermi gas. The location of the kink (Fermi surface when A = 0) 
is determined by Luttinger theorem. 

For a Hubbard model with a strong on-site attraction, the ground state 
is the Bose condensate of local pairs. This corresponds to a p independent 
6p and the corresponding Fermi distribution function is a plateau within the 
Brillouin zone. 

In case of the boson-fermion model, the distribution function corresponds 
to an different interpolation from that in [13] between the former two cases. 
The np profile is the superposition of the BCS case with the Fermi level 
retreated towards the bottom of the band and with a plateau outside the 
Fermi sea. 

The Fermi distribution function for the three different cases is shown 
schematically in Fig. 2. The profile in Fig. 2c may be generated with 
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Figure 2 The schematic Fermi distribution function at T = for BCS model (a), local 
pair model (b) and boson-fermion model (c), where the gap smearing effect is not shown. 

a Hubbard like model with an on-site attraction and a nearest neighbor 
repulsion [14]. 

5.2 The pole trajectory underneath the physical sheet at T > Tq 

In the BCS model, the pairing instability would be triggered at the ab- 
sence of a long range order by a pair of complex conjugate poles of the two 
fermion scattering amplitude below Tq on the energy .E-plane. For T > Tc, 
they slip through the cut at E = lap to the unphysical sheets below and 
leave the real axis vertically with an increasing T. Their effect becomes in- 
significant outside the critical region. On the other hand, for the local pair 
system the bound state pole stays always on the real E axis for T > Tc- 
The pole trajectory for the boson-fermion model (2.1) (for which the pole 
of the scattering amplitude of two fermions coincides with the pole of the 
boson propagator ) has been calculated to the one-loop order. Without a 
Bose condensate, there would be a pair of complex conjugate poles trigger- 
ing the instability for T < Tc- As T being increased away from T = Tc, 
they emigrant to the unphysical sheet but leave the real axis obliquely with 
a tunable slope proportional to g 2 (Fig. 3). For small g, they stay close to 
the real E for a wide range of temperature above Tc and could make signif- 
icant contributions to transport processes. Starting with a purely electronic 
Hamiltonian, such resonance boson poles can be implemented for a few body 
system [14]. But to reproduce the pole trajectory of the boson-fermion model 
for a many-body system remains an open problem. 

I am grateful to Professor R. Micnas and the organizing committee for 
inviting me to this prestigious conference. I am also benefited from many 
discussions with Professors R. Micnas and S. Robaszkiewicz. This work is 
supported in part by National Science Council of ROC under Grant NSC- 
CTS-981003 and by U. S. Department of Energy under Grant DE-FG02- 
91ER40651, Task B. 
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Figure 3 The pole trajectory on the complex energy plane of the two electron (hole) 
scattering amplitude for T > Tq- The dashed line is for the boson- fermion model and 
the dotted line is for BCS model. The arrows point to the direction corresponding to an 
increasing T The heavy solid line represents the branch cut and the intersection point of all 
trajectories corresponds to E = 2cf- 
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